MEASUREMENT OF THE TOP QUARK MASS 



STEVEN R. BLUSK 
Department of Physics and Astronomy, 
University of Rochester, 
Rochester, NY, 14627 

for the 

CDF and DO Collaborations 




The first evidencea and subsequent discoverjoCi of the top quark was reported nearly 4 years 
ago. Since then, CDF and DO have analyzed their full Run 1 data samples, and analysis 
techniques have been refined to make optimal use of the information. In this paper, we report 
on the most recent measurements of the top quark mass, performed by the CDF and DO 
collaborations at the Fermilab Tevatron. The CDF collaboration has performed measurements 
of the top quark mass in three decay channels from which the top quark mass is measured to 
be 175.6±6.8 GeV/c^. The DO collaboration combines measurements from two decay channels 
to obtain a top quark mass of 172.1 ±7.1 GeV/c^. Combining the measurements from the two 
experiments, assuming a 2 GeV/c^ correlated systematic uncertainty, the measurement of the 
top quark mass at the Tevatron is 173.9 ±5.2 GeV/c^. This report presents the measurements 
of the top quark mass from each of the decay channels which contribute to this measurement. 



1 Introduction 



The top quark is defined as the I3 = +1/2 member of a weak SU(2) isodoublet that also contains 
the b quark. In pp collisions, top quarks are produced primarily in ti pairs and are expected to 
decay through the electroweak interaction to a W-boson and b-quark. In the standard model, 
the branching fraction for t Wb is expected to be nearly 100%. The decay width§ is calculated 
to be 1.6 — 1.7 GeV for masses between 150 - 180 GeV/c^. The top quark mass is sufficiently 
large that top-flavored hadrons are not expected to formQ. 

The top quark mass, Mf^p , is an important parameter in calculations of electroweak processes 
since its mass is approximately 35 times larger than that of the next heaviest fermion. Like other 
fermions, its mass is taken to be an unpredicted parameter of the standard modeli. Calculations 
of electroweak radiative corrections relate the top quark and the W-boson {M]y) masses to that 
of the Higgs boson. Precision measurements of Mtop and Myy can therefore aid in searches for 
the Higgs boson. 



2 Overview of Top Quark Mass Measurement 

In ^/s = 1.8 TeV pp collisions, top quarks are expected to be produced in pairs mainly through 
the subprocess qq — > ti. At the Tevatron, both CDF and DO select samples of events which are 
consistent with containing a tt pair. Since each W-boson may decay leptonically or hadronically, 
this leads to three categories of events. Events which have 2 high Pt isolated leptons (both W- 
bosons decay leptonically), missing transverse energy (HSt )^ and two or more jets (from the 
b-quarks) are classified as Dilepton events. Events which have one isolated high Pt lepton, 
Eix , and 4 or more jets are referred to as l+jets events. Events containing 6 or more jets and 
no high Pt isolated leptons are classified as All-hadronic events. Neither CDF nor DO try to 
reconstruct the top mass from candidate ti ^ t + X events. Ignoring final states containing r's, 
the branching ratio for ti into these three decay channels is 5%, 30%, and 44% respectively. 

Both CDF and DO have measured the top quark mass in the Dilepton and l+jets channels. 
CDF has also made a measurement of the top quark mass in the All-Hadronic channel. A 
general procedure is used for extracting the top quark mass. For each decay channel, we choose 
a variable (or a set of variables) which can be measured in the data sample, and is highly 
correlated with Mtop ■ Examples of such variables include event-by-event top mass, jet Et, 
Fix ) total scalar Et, etc. The distribution of a given variable is modelled using Monte Carlo 
(MC) simulations of ti signal and background. The ti signal distributions are modelled using 
the HERWIci MC simulation and are produced for a wide range of input top quark masses, 
ranging from ~120 to 220 GeV/c^. Background distributions are derived from a combination 
of Monte Carlo and/or data, depending on the channel. The modelled signal and background 
distributions are generically referred to as templates. The top quark mass is extracted from a 
maximum likelihood procedure which compares the observed data distribution(s) to the signal 
and background templates. Assumptions made in the modelling of signal and background are 
included as sources of systematic uncertainty. 



3 Top Quark Mass Measurement in the Dilepton Topology 

Both CDii and D(0 have reported measurements of the top quark mass in the Dilepton decay 
topology. Because these events are presumed to have two neutrinos, the kinematics of the 
ti system (including the top quark mass) cannot be uniquely determined. A measurement of 
the top quark mass can however be made by comparing the measured observables with the 
expectation from ti and background. Each experiment employs two techniques to measure the 
top quark mass from their respective samples. 



3.1 CDF Dilepton Mass Measurements 



The first Dilepton mass measurement uses the energies of the 2 highest Ex jets in the can- 
didate events as a discriminator for the top quark mass measurement. Prom a sample of 8 
events, with an expected background of l.libO.3 events, the analysis yields a top quark mass of 
159ib23(stat)ibll(syst) GeV/c^. A second analysis exploits the fact that, in the W-boson rest 
frame of the t — > Wb decay, the top quark mass can be related to the invariant mass of the b-jet 
and lepton {Mi},) to which it decays. Among the two possible pairings of the 2 leptons with the 2 
highest Et jets, the pair with the smaller sum of invariant masses is chosen. From this analysis, 
the top quark mass is measured to be 163±20(stat)±6(syst) GeV/c^ The jet energy scale, and 
signal and background modelling dominate the systematic uncertainty. The two measurements 
are combined taking into account the correlations in uncertainties to obtain a top quark mass 
of 161±17(stat)±10(syst) GeV/c^. 

3.2 DO Dilepton Mass Measurements 

Both techniques which DO employs uses MC simulations to calculate a probability that the ob- 
served kinematics of the event are consistent with top of an assumed mass. This probability 
is mapped out as a function of assumed top mass in the range from 80 to 280 GeV/c^. The 
probability in the first technique is based on: (a) the probability for the charged leptons to have 
come from tt decay at the assumed mass, and (b) the probability for the event to have been 
produced by valence quarks with momentum fractions f{x) and /(x). The second technique 
calculates a probability based on the difference between the measured and calculated in 
the event. The calculation of the EIt is facilitated by assuming values for the pseudorapidity 
of the two neutrinos (r/i and r/2), and then integrating over the rji, 772 phase space. For both 
techniques, detector resolution effects are included as well as trying both lepton- (b-quark) com- 
binations. Based on the data sample of 6 events and an expected background of 1.42 events, 
the mass is measured to be 168.1±12.4(stat)±3.6(syst) GeV/c^ using the first technique, and 
169.9±14.8(stat)±3.6(syst) GeV/c^ using the second technique. The dominant systematics are 
the same as those mentioned for the CDF dilepton analyses. The measurements are combined 
to obtain a top quark mass of 168.4ibl2.3(stat)±3.6(syst) GeV/c^. 

4 Top Quark Mass Measurement in the l-(-jets Topology 

Both CDI§ and Dtltl have measured the top quark mass in the 1+jets topology. The 1+jets 
decay channel currently allows for the most accurate determination of the top quark mass. The 
measurement benefits from a relatively large branching ratio and the ability to fully reconstruct 
the top mass on an event-by-event basis. The momentum of all final state particles are assumed 
to be measured, with the exception of the unobserved neutrino, v. The transverse momentum 
components of the neutrino however can be inferred from the in the event, leaving only 
its Z-component unmeasured. Three kinematic constraints are implied by the tt — > IvqqhbX 
decay hypothesis. Namely, we require that Miy = M]y, Mjj = M\y, and Mi^j, = Mjjh, where 
b = b OT b jet, j=non-6 jet, Mt^=W-boson mass, and -Mj,a,(2,) are invariant masses formed 
from the reconstructed 4-vectors of the indicated objects. Because ti events are expected to 
contain 2 b-quarks, as compared to only a few percent for background events, identifying b-jets 
improves the signal to background ratio (S/B). B-jets are tagged by either reconstructing the 
secondary decay vertices of B-hadrons (SVX tagged jet), or by identifying charged tracks which 
are consistent with coming from semileptonic b-hadron decay (SLT tagged jet). Each event is fit 
to the ti — > luqqbbX topology, requiring b-tagged jets to be assigned to b-quarks. All allowable 
permutations are tried, and the top mass corresponding to the solution with the lowest < 10 
is taken as the reconstructed top mass (Mrec) for a given event. The ensemble of reconstructed 



masses (one per event) are fit using a likelihood procedure to extract the best estimate of the 
top quark mass. 

4.1 CDF l+jets Mass Measurement 

The initial sample selection requires 1 isolated lepton with Pt > 20 GeV/c, |r/| < 1, El^ > 
20 GeV, > 3 jets with Et > 15 GeV, |r/| < 2.0, and a fourth jet with Et > 8 GeV and |ry| < 2.4. 
After applying these selection criteria and the cut on the kinematic fit, the sample consists of 
151 events, of which ~35% are estimated to be from tt. Motivated by the different S/B and top 
mass resolution for events with SVX, SLT, or no tags, MC simulations were used to determine 
an optimal wav to partition the sample as to obtain the best precision on the top quark mass 
measurements. Studies showed that an optimal partitioning consists of subdividing the events 
into the four statistically independent categories shown in Table!. The table shows the numbers 
of events in each subsample, the expected background fractioicj Xb, and the fitted top mass. 



Data Sample 


T^Events 




Top Mass (GeV/c^) 


SVX Double (2 SVX tagged jets) 
SVX Single (1 SVX tagged jet) 
SLT tagged (>1 SLT tagged jets, no SVX tags) 
No Tags (>4 jets with Et > 15 GeV, < 2.0) 


5 

15 
14 
42 




170.1 ±9.3 
178.0 ±7.9 

180.8 ±9.0 


Combined 


76 




175.9 ±4.8 



Table 1: Subsamples used in the top quark mass analysis, expected background fractions, and the measured top 
quark mass in each sample. The last line shows the results from the combined likelihood fit to all four subsamples. 



Since the four subsamples are statistically independent, the joint likelihood for the 76 events 
is given by the product of the four subsample likelihoods. Each subsample likelihood uses 
parameterized signal and background probability densities, along with a constraint on the the 
background fraction (see Table 1), to evaluate the likelihood for observing the measured set of 
masses as a function of the top mass. The results of the likelihood fit for each subsample and the 
combined result are shown in Table 1 . The results of the combined fit for the 76 data events and 
are shown in Fig. 1. Including the systematic uncertainties presented in Table 3, the top quark 
mass is measured to be 175.9±4.8(stat.)±4.9(syst.) GeV/c^ for events in the l±jets topology. 

4-2 DO l+jets Mass Measurement 

The event selection for the DO analysis is similar to CDF. The main differences are a larger 
\r]\ cut for the charged lepton (from the W-boson decay) and all four jets are required to have 
Et > 15 GeV. After the event selection criteria and the kinematic cut, 77 events remain 
of which 5 are SLT tagged (DO did not have a silicon vertex detector in Run 1). DO extends 
the likelihood procedure by including an additional variable D which provides discrimination 
between signal and background events. Two versions of the discriminant are tried. The first 
is a low-bias version Dn,, and the second, D^j^f^ was the output from a neural network trained 
on ti signal and background events. Both discriminants provide similar levels of discrimination 
between signal and background events. Two-dimensional signal templates of D vs M^ec are 
generated for many input values for the top mass and a single template for the background. 
For each top mass value, a likelihood is minimized with respect to the number of signal and 
background events. The corresponding negative log-likelihood values are fit to a quadratic 
function of top mass, and the measured top mass is taken to be where the function is a minimum. 
The fits to the data are presented in Table 2 and shown in Fig. 2 For illustrative purposes, the 
data sample is separated into (a) signal depleted and (b) signal enhanced subsamples. Also 



shown are the hkehhood fits for the two methods. Including the systematic uncertainties shown 
in Table 3, the top quark mass is measured to be 173.3±5.6(stat.)±5.5(syst.) GeV/c^. 



Method 


Mtop (GeV/c^) 


Fitted Ns 


Fitted Nb 


LB 
NN 
LB+NN 


174.0±5.6 
17L3±6.0 
173.3±5.6 


24±8 
29±9 


53±10 
48±10 



Table 2: Fit results using the low-bias (LB) likelihood, then neural-network (NN) method, and the combined 
result. Shown are the fitted top quark masses, and the fitted number of signal {Ng) and background {Nb) events. 



Source 


CDF 


DO 


Energy Scale 


4.4 


4.0 


Signal Model 


1.8 


1.9 


Background 


1.3 


2.5 


Noise/M.L 




1.3 


b-tagging bias 


0.4 




Likelihood Method 




1.3 


Total 


4.9 


5.5 



Table 3: Main systematic uncertainties (GeV/c^) in the l-|-jets top quark mass measurement for each experiment. 
(M.I.=Multiple Interactions; PDF's = Parton Distribution Functions) 

5 CDF AU-Hadronic Top Mass Measurement 

This measuremenlEi requires 6 or more jets with Et > 15 GeV, and \rj\ < 2, and at least 
one b-tagged jet. Additional cuts, based on MC simulations, were employed to improve the 
S/B ratio. An event-by-event top mass is evaluated using the measured jet energies and the 
kinematic constraints mentioned previously. All thirty combinations are tried, and the top mass 
having the lowest < 10 is taken as the top mass for the event. The data sample consists of 
136 events, of which 108+9 events are expected to come from background sources. A likelihood 
procedure (with a gaussian background constraint) was used to extract a measured top mass of 
186±10(stat.)±12(syst.) GeV/c^. The systematic uncertainties are dominated by the jet energy 
scale signal modelling, and the fitting method. 

6 Combining the top mass measurements 

The measurement in each decay topology may be combined to improve the resolution on the top 
quark mass measurement. Each experiment combines their measurements into a single measure- 
ment, taking into account correlations in uncertainties between the measurements. Combining 
the three CDF measurements, the top quark mass is measured to be 175.6ib6.8 GeV/c^. The two 
DO measurements yield a top quark mass of 172.1ib7.1 GeV/c^. These two mass measurements 
are combined assuming a 2 GeV/c^ correlated systematic uncertainty. The top quark mass at 
the Tevatron is therefore measured to be 173.9±5.2 GeV/c^. The values of Mw versus Mtop for 
several (standard model) Higgs masses are shown in Fig. 3 The data tend to favor a mass below 
SS500 GeV/c^, but clearly more data is needed to reduce the uncertainties. In Run 2, we expect 
to measure Mtop to within 1-2 GeV/c^ and Mw to within 40 MeV/c^. This will clearly narrow 
the search window for future Higgs' searches. 



I would like to thank my collaborators at CDF for their assistance in preparing this talk. 
I would also like to thank the DO top group, particularly Scott Snyder and Erich Varnes for 
answering questions regarding their top mass analyses. 
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Figure 1: The results of applying the likelihood procedure to the combined subsamples. The figure shows the 
data (points), fitted background (light shaded region), and fitted signal plus background (darker shaded region). 
The inset shows the shape of — log L versus the generated top mass from which we extract the best estimate of 

the top mass and the statistical uncertainty. 
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Figure 2: The fit results from tfie LB and NN likeliliood procedures applied to the 77 event data sample, (a) 
Events which have Dw > 0.43 and (b) Dib < 0.43. The histogram shows the data, and the points are the MC fit 
results, (c) shows the fit to —logC vs top mass for the two methods. 




Figure 3: Mw versus Mtop for various choices of the Higgs mass. Shown are the CDF and DO averages, as weU 

as the world average. 



